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Introduction 

Applications of direction finding 

While direction finding for navigation purposes (referred to 
as cooperative direction finding) is becoming less impor¬ 
tant due to the availability of satellite navigation systems, 
there is a growing requirement for determining the loca¬ 
tion of emitters as the mobility of communications equip¬ 
ment increases: 

i In radiomonitoring in line with ITU guidelines 

■ Searching for sources of interference 

■ Localization of non-authorized transmitters 
i In security services 

■ Reconnaissance of radiocommunications of criminal 
organizations 

i In military intelligence [1] 

■ Detecting activities of potential enemies 

■ Gaining information on enemy's communications order 
of battle 

i In intelligent communications systems 

■ Space division multiple access (SDMA) requiring 
knowledge of the direction of incident waves [2] 


i In research 

■ Radioastronomy 

■ Earth remote sensing 

Another reason for the importance of direction finding lies 
in the fact that spread-spectrum techniques are increas¬ 
ingly used for wireless communications. This means that 
the spectral components can only be allocated to a spe¬ 
cific emitter if the direction is known. Direction finding is 
therefore an indispensable first step in radiodetection, par¬ 
ticularly since reading the contents of such emissions is 
usually impossible. 

The localization of emitters is often done in a multistage 
process. Direction finders distributed across a country 
allow an emitter to be located to within a few kilometers 
(typ. 1 % to 3% of the DF distances) by means of triangu¬ 
lation. The emitter location can more precisely be deter¬ 
mined with the aid of direction finders installed in vehicles. 
Portable direction finders moreover allow searching within 
the last 100 m, for instance in buildings. 

Historical development 

The DF technique has existed for as long as electromag¬ 
netic waves have been known. It was Heinrich Hertz who 
in 1888 found out about the directivity of antennas when 
conducting experiments in the decimetric wave range. A 
specific application of this for determining the direction of 
incidence of electromagnetic waves was proposed in 1906 
in a patent obtained by Scheller on a homing DF method. 
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The initial DF units were polarization direction finders. 

They consisted of a rotatable electric or magnetic dipole 
whose axis was brought to coincidence with the direc¬ 
tion of the electric or magnetic field. From the direction of 
polarization, the direction of incidence was then deduced. 
The rotating-loop direction finder is one of the best known 
direction finders of this type. In 1907, Bellini and Tosi 
discovered the DF principle that was named after them: a 
combination of two crossed directional antennas (e.g. loop 
antennas) with a moving-coil goniometer for determining 
the direction [1]. Despite this invention, rotating-loop 
direction finders were often used in the First World War 
(Fig. 1). 

The invention made by Adcock meant a great leap forward 
in the improvement of the DF accuracy with respect to 
skywaves in the shortwave range. The pharmacist by pro¬ 
fession realized in 1917 that with the aid of vertical linear 
antennas (rod antennas or dipoles) directional patterns can 
be generated that correspond to those of loop antennas 
but do not pick up any interfering horizontally polarized 
field components (G. Eckard proved in 1972 that this is not 
true in all cases [3]). It was not until 1931 that Adcock an¬ 
tennas were first employed in Great Britain and Germany. 

In 1925/26, Sir Watson-Watt took the step from the me¬ 
chanically moved goniometer direction finder to the elec¬ 
tronic visual direction finder. As from 1943, British naval 
vessels were equipped with crossed loops and three- 
channel Watson-Watt direction finders for the shortwave 
range ("huff-duff" for detecting German submarines). 

As from 1931, camouflaged direction finders were avail¬ 
able for use in vehicles and as portable direction finders 
for detecting spies. 


The first shortwave direction finder operating on the 
Doppler principle was built in 1941. The rapid progress 
in the development of radar in Great Britain made it 
necessary to cover higher frequency ranges: In 1943, the 
first direction finders for "radar observation" at around 
3000 MHz were delivered [1 ]. 

As from 1943, wide-aperture circular-array direction find¬ 
ers (Wullenweber) were built for use as remote direction 
finders. Since the 1950s, airports all over the world have 
been equipped with VHF/UHF Doppler direction finding 
systems for air traffic control. 

In the early 1970s, digital technology made its way into 
direction finding and radiolocation; digital bearing genera¬ 
tion and digital remote control are major outcomes of this 
development. 

Since 1980, digital signal processing has been increasingly 
used in direction finding. It permits the implementation of 
the interferometer direction finder and initial approaches 
toward the implementation of multiwave direction finders 
(super-resolution). The first theoretical considerations were 
made much earlier, e.g. in [4]. 

Another important impetus for further development 
came from the requirement for the direction finding of 
frequency-agile emissions such as frequency-hopping and 
spread-spectrum signals. The main result of this was the 
broadband direction finder, which is able to simultaneous¬ 
ly perform searching and direction finding based on digital 
filter banks (usually with the aid of fast Fourier transforms 
(FFTs)) [5]. 



Fig. 1: Mobile 
rotating-loop direction 
finder for military use 
(about 1918) [1], 
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Tasks of direction finding 

The task of a radio direction finder is to estimate the direc¬ 
tion of an emitter by measuring and evaluating electro¬ 
magnetic field parameters. 

Usually, the azimuth a is sufficient to determine the di¬ 
rection; the measurement of elevation e is of interest for 
emitters installed on flying platforms and especially for the 
direction finding of shortwave signals (Fig. 2). 

Only in the case of undisturbed wave propagation is the 
direction of the emitter identical with the direction of inci¬ 
dence of the radio waves. Usually, there is a large number 
of partial waves arriving from different directions and mak¬ 
ing up a more or less scattered field. The direction finder 
takes from this wavefront spatial and temporal samples 
and supplies in the ideal case the estimated values a and 
s as the most probable direction of the emitter observed. 

Bearings can be taken using the following reference direc¬ 
tions (Fig. 3) (see also DIN EN3312 [6]): 
i Geographical north (true north) > true radio bearing 
i Magnetic north 

i Vehicle axis > relative or direct radio bearing 


DF principles 

Generation and characterization of electromagnetic 
waves 

Electromagnetic waves are caused by charging and dis¬ 
charging processes on electrical conductors that can be 
represented in the form of alternating currents [7], [8]. 

The first assumption is based on the undisturbed propaga 
tion of a harmonic wave of wavelength A. At a sufficiently 
large distance, the radial field components are largely 
decayed so that, limited to a small area, the wave can be 
considered to be plane: Electric and magnetic field com¬ 
ponents are orthogonal to and in phase with one another 
and are perpendicular to the direction of propagation, 
which is defined by the radiation density vector (Poynting 
vector) S 



Z 

^“0 

where E = RMS value of electric field strength 

Z 0 = characteristic impedance of free space; Z 0 = 1207iQ 

or by the wave number vector k (Fig. 4). 


Fig. 2: Definition of emitter direction 



Fig. 3: Reference directions 


North 
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Overview of main DF principles 

Direction finding relies on the basic characteristics of elec¬ 
tromagnetic waves, which are: 

i Transversality, i.e. field vectors are perpendicular to the 
direction of propagation 

i Orthogonality of phase surfaces and direction of 
propagation 

Every DF process essentially employs one of the following 
methods (Table 1): 

i Method A: measures the direction of electric and/or 
magnetic field vectors > polarization direction finders 
i Method B: measures the orientation of surfaces of equal 
phase (or lines of equal phase if the elevation is not of 
interest) > phase direction finder 

Polarization direction finders are implemented by means of 
dipole and loop antennas. The classic rotating-loop direc¬ 
tion finder also belongs to this category (rotation of loop to 
minimum of received signal > direction of incident wave 
perpendicular to loop). Today, polarization direction find¬ 
ers are used in applications where there is sufficient space 
only for small antennas, e.g. in vehicles and on board 
ships for direction finding in the HF band. Evaluation is 
usually based on the Watson-Watt method (see section 
"Classic DF methods"). To obtain unambiguous DF results, 
however, the phase must be evaluated in addition: 


Fig. 4: Propagation of space waves 



Unambiguous DF results require evaluation of the phase 
information. 

Phase direction finders obtain the direction information 
(bearing information) from the spatial orientation of lines 
or surfaces of equal phase. There are two basic methods: 

i Direction finding based on directional patterns 
With this method, partial waves are coupled out at vari¬ 
ous points of the antenna system and combined at one 
point to form a sum signal. The maximum of the sum 
signal occurs at the antenna angle at which the phase 
differences between the partial waves are at a minimum. 
The sum signal is thus always orthogonal to the phase 
surfaces of the incident wave (maximum-signal direction 
finding). For minimum-signal direction finding, the partial 
waves are combined so that the phase differences in the 
direction of the incident wave become maximal and there 
is a distinct minimum of the received signal 
i Direction finding by aperture sampling 
With this method, samples are taken at various points of 
the field and applied to evaluation circuits sequentially or 
in parallel. These circuits determine the bearing by link¬ 
ing the samples together, which is today mostly done by 
mathematical operations. 

Typical examples are interferometers and Doppler direction 
finders. 

The DF methods mentioned so far are suitable only to a 
limited extent for determining the directions of incidence 
of several waves that overlap in the frequency domain. 


Table 1: DF principles. 


Wave characteristic 

Transversality 

Phase surfaces 

Direction of propagation 

DF principle 

Polarization direction 
finder 

Phase direction finder 

Examples 


Direction finding with directional 
patterns 

Aperture sampling 



Conversion phase > amplitude 

Direct evaluation 

Sensor array processing 


1 Rotating loop 

1 Directional antenna 

Maximum-signal direction finder 
Minimum-signal direction finder 

1 Interferometer 

1 Correlation direction 
finder 


1 Dipole 

1 Adcock with Watson-Watt 
evaluation 

1 Rotating field direction 
finder 

1 Adaptive beam former 


1 Loaded loop 


1 Doppler direction 
finder 

1 MUSIC 


1 Crossed loop with 
Watson-Watt evaluation 



1 ESPRIT 
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With the progress made in digital signal processing, the 
methods known from the theory of spectral estimation 
have been applied to the analysis of wavefront and further 
developed. The term "sensor array processing" describes 
the technique of gaining information about the parameters 
of incident waves from the signals derived from the ele¬ 
ments, or sensors, of a sensor array (antenna array in di¬ 
rection finding, hydrophone array for sonar). 

There are basically three different methods: 
i Beamforming methods, e.g. correlation direction finder, 
spatial Fourier analysis, adaptive antenna 
i Maximum likelihood method as the most general 
model-based method 
i Subspace methods, e.g. MUSIC, ESPRIT 

Main requirements on DF systems: 

i High accuracy 
i High sensitivity 
i Sufficient large-signal immunity 
i Immunity to field distortion caused by multipath 
propagation 

i Immunity to polarization errors 
i Determination of elevation in shortwave range 
i Stable response in case of non-coherent co-channel 
interferers 

i Short minimum required signal duration 
i Scanning direction finders: high scanning speed, and 
high probability of intercept (POI) 



- level 

- band occupancy 


Components of a DF system 

A DF system (Fig. 5) consists of the following compo¬ 
nents: 

i Antenna system 
i DF converter 
i Evaluation unit 
i Display unit 

Depending on the configuration, systems for determining 
the direction finder's own coordinates/orientation (GPS, 
compass), remote-control units (LAN, WAN), antenna con¬ 
trol units, etc, can be added. 

The achievable DF speed mainly depends on the number 
H of receive sections, as this parameter determines the 
number of antenna outputs that can be measured in paral¬ 
lel. 

To achieve maximum speed, it must be possible to gen¬ 
erate a bearing in a single time step, i.e. from one set of 
samples (monopulse direction finding). For unambigu¬ 
ous direction finding over the total azimuth range, at least 
three antenna outputs are required. If there are also three 
receive sections, multiplexing of the measurement channel 
is not necessary. 

Typical examples of monopulse DF antennas: 
i Multimode antenna for amplitude comparison direction 
finders, e.g. Adcock antenna 
i Interferometer and rotating-field (phase) direction finder 

For high DF accuracy (e.g. 1°) and large bandwidth (e.g. 

1 MHz to 30 MHz or 20 MHz to 1000 MHz), five to nine 
aperture samples are usually required. Since monopulse 
solutions would then be very complex, one fixed and two 
sequentially switched receive sections are frequently used. 

The DF converter converts the carrier-frequency antenna 
signals to a fixed IF. Since this conversion must be per¬ 
formed with equal phase and amplitude in all receive 
sections, the use of a common synthesizer is indispens¬ 
able. Moreover, with most multipath direction finders, the 
receive sections are calibrated in order to ensure equal 
amplitude and phase. This is done with the aid of a test 
generator at defined intervals and prior to the actual DF 
operation. 

The evaluation unit determines the bearing from the ampli¬ 
tudes and/or phases of the IF signal. 
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Classic DF methods 

Using directional antennas 

Evaluating the receive voltage of a mechanically rotated 
directional antenna with reference to the direction is the 
simplest way of direction finding. With this method, the 
bearing is derived from the characteristic of the receive 
voltage as a function of the antenna rotation angle: When 
a wave arrives, the receive voltage yields the directional 
pattern of the antenna. The pattern position relative to the 
antenna rotation angle is the measured bearing [1]. 

This type of direction finder is a phase direction finder 
since the directivity of the receive antenna is achieved by 
superimposing partial waves whose phase differences 
depend on the angle of incidence. In the simplest case, 
the antenna is rotated and the bearing determined by the 
operator. The antenna is rotated until the receiver output 
voltage assumes an extreme value. The antenna direction 
thus found is read from a scale, and the bearing is deter¬ 
mined from it. If the directional antenna (with maximum 
or minimum pattern) is permanently rotated with the aid 
of a motor, and the receive voltage is displayed graphically 
as a function of the angle of rotation, a rotating direction 
finder is obtained (Fig. 6). Using suitable automatic evalua¬ 
tion of the receive voltage characteristic, e.g. by means of 
a maximum detector, a fully automatic direction finder is 
obtained. 



The following benefits are common to all variations of this 
DF method: 

i High sensitivity due to the directivity of the antenna 
i Simple and inexpensive implementation (only one 
receiver required (single-channel principle)) 
i Resolution of multiwavefronts is possible (prerequisite: 
different angles of incidence and high-directivity antenna 
system) 

i Same antenna can be used for direction finding and 
monitoring 

The drawbacks of this method result from the restricted 
angular detection range, which is due to the directivity 
of the antenna, and the antenna's limited rotating speed, 
which is mainly due to the use of a mechanical rotator: 
i Probability of intercept is reciprocal to the directivity 
i Method fails in case of short-duration signals, i.e. with 
signal dwell times that are short compared to a scanning 
cycle of the antenna 

Despite these drawbacks, DF methods using mechanically 
rotated directional antennas are still in use today since 
achieving the described advantages with other methods 
involves in part considerably higher cost and effort. In the 
microwave range, in particular, the mechanical DF method 
is often the only justifiable compromise between gain, low 
noise and expenditure. 

If, in addition to a directional pattern with a maximum in 
the direction of the incident wave, a directional pattern 
with a minimum is used, a monopulse direction finder is 
obtained that even with a slowly rotating or fixed antenna 
furnishes bearings as long as waves arrive within the main 
receiving direction of the antenna. Fig. 7 shows an imple¬ 
mentation of this using log-periodic dipole antennas that 
are connected by means of a 0/180° hybrid. The directional 
patterns shown below are thus obtained. 


Fig. 7: DF using sum-difference method 




Typical implementation on the left, directional patterns for sum (I) and 
difference (A) outputs on the right. 
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The quotient of the difference and the sum signal yields 
a dimensionless, time-independent function, i.e. the DF 
function: 


PF(a) = 


Ua(«) 

U z (a) 


After forming the quotient of the two test voltages, the DF 
function immediately furnishes the bearing a. 


If there is a phase difference 5 between the two voltages 
U x and U y , which may be due to ambient interference (e.g. 
reflections), the displayed figure is an ellipse. The position 
of the main axis yields the bearing a , which is calculated 
from the two voltages by means of the equation below [9], 


r 


a = Re 


arctan 


A 


v 


U 


1 2U, 

= —arctan- 


U. 


cos 5 


y J 


U, 


U 


Watson-Watt principle 

If the amplified and filtered signals of a receiving antenna 
with outputs for a sine-shaped and a cosine-shaped direc¬ 
tional pattern are applied to the x and y deflection plates 
of a cathode-ray tube (CRT), a line Lissajous figure is ob¬ 
tained in the ideal case, whose inclination a corresponds 
to the wave angle but exhibits an ambiguity of 180°. The 
indicated angle is obtained from the ratio of the two sig¬ 
nals as follows: 

+ u x 

a = arctan— 

U y 

An unambiguous bearing indication is obtained (Fig. 8) if 
a blanking signal is additionally used in this DF method, 
which was first implemented by Watson-Watt in 1926. The 
blanking signal is derived from an omnidirectional receiv¬ 
ing antenna with an unambiguous phase relationship. 


The principal advantage of this method is that the bearing 
is indicated without delay, which means that it is capable 
of monopulse direction finding over the entire azimuth 
range. 

Suitable antennas (Fig. 9) with sine-shaped or cosine¬ 
shaped directional patterns are in particular the following: 
i Loop antennas (or ferrite antennas) 
i Adcock antennas (monopole or dipole arrays) 

Crossed-loop antennas with Watson-Watt evaluation are 
mainly suitable for mobile applications due to their com¬ 
pact size. They feature the following benefits and draw¬ 
backs: 

Benefits 

i Extremely short signal duration is sufficient 
i Implementation is simple 
i Minimum space is required 


Fig. 8: Watson-Watt direction finder with crossed-loop antenna 



Blanking 
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Drawbacks 

i Small-aperture system (D/A < 0.2) causing errors in case 
of multipath propagation 

i Large DF errors when receiving sky waves with steep 
elevation angles 

Adcock antennas feature the following advantages over 
crossed-loop antennas: 

i Improved error tolerance for skywave reception 
i Wider apertures can be implemented to reduce errors 
in case of multipath reception (e.g. D/A < 1 for 8-fold 
Adcock) 

Modern direction finders no longer display the IF voltages 
of antenna signals on a CRT but digitally process the sig¬ 
nals after converting them into a relatively wide IF band 
(Fig. 10). 



Configuration of a modern direction finder operating on the 
Watson-Watt principle. 


Selection is mainly effected by means of digital filters; 
bearings are calculated numerically, e.g. using the last 
equation above, and displayed on a computer (worksta¬ 
tion, PC) with a graphical user interface (GUI). 

A number of disadvantages of analog direction finders are 
thus avoided, yielding the following effects: 
i Synchronous operation of channels also on filter edges 
i Simple procedure of taking into account correction 
values for antenna networks, cables, etc 
i No temperature drift in digital section 
i Bearings available as numeric values for further 
evaluation, especially for easy transmission to remote 
evaluation stations 


Doppler direction finder 

If an antenna element rotates on a circle with radius R, the 
received signal with frequency co 0 is frequency-modulated 
with the rotational frequency co r of the antenna due to the 
Doppler effect. If the antenna element moves toward the 
radiation source, the receive frequency increases; if the 
antenna element moves away from the radiation source, 
the receive frequency decreases. 


From the instantaneous amplitude 


f 


2ttR 


A 


u(t)= acos((j)(t))= a cos oo 0 t +-cos(co r t-a)+cp 

hr 


V 


J 


the instantaneous frequency is derived by differentiation of 
the phase 

d<|>(t) 27iR . , x 

co(tj = —— = co 0 -co r sin(co r t-ocj 


dt 


A, 


After filtering out the DC component go 0 , the demodulated 
Doppler signal is obtained as 


27iR . , x 

S D =——co r sin(co r t-ocj 


A 


0 


Fig. 9: Various antenna configurations 



For generating sine-shaped or cosine-shaped directional patterns: crossed-loop antenna, ferrite antenna, H Adcock, U Adcock (from left). 
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The phase of the demodulated signal is compared with a 
reference voltage of equal center frequency derived from 
the antenna rotation 

S r = — sin cD r t 

which yields the bearing a [1]. 

Since rotating an antenna element mechanically is nei¬ 
ther practically possible nor desirable, several elements 
(dipoles, monopoles, crossed loops) are arranged on a cir¬ 
cle (Fig. 11) and electronically sampled by means of diode 
switches (cyclic scanning). 

To obtain unambiguous DF results, the spacing between 
the individual antenna elements must be less than half the 
operating wavelength; in practice, a spacing of about one 
third of the minimum operating wavelength is selected. 

If this rule is adhered to, Doppler DF antennas of any size 
are possible, i.e. wide-aperture systems featuring the fol¬ 
lowing can easily be implemented: 
i High immunity to multipath reception 
i High sensitivity 

A drawback of the Doppler method is the amount of time 
it requires. At least one antenna scanning cycle is required 
in order to obtain a bearing. At a typical rotational fre¬ 
quency of 170 Hz in the VHF/UHF range, one cycle takes 
approx. 6 ms. 

Interferometer 

The interferometer direction finder was first used in radio 
astronomy [10]. The objective there was to increase the 
resolution power and the sensitivity of the DF system by 
superimposing the signals of only a few antenna elements 
that were however spaced many wavelengths apart 
(Fig. 12). 


The ambiguous interference pattern is weighted with the 
directional pattern of the antenna element, which yields an 
unambiguous bearing (see Fig.12, diagrams on the right). 

Using the minimum of three omnidirectional receiving el¬ 
ements, unambiguous determination of the azimuth and 
elevation is possible only if the spacing a between the an¬ 
tennas is no greater than half the wavelength. With 0 1f cD 2 , 
O 3 as the phases measured at the outputs of the antenna 
elements, the azimuth is calculated as 

a = arctan——— 

The elevation angle is obtained as 

(q 2 -O i ) 2 + ( ( D3-O i ) 2 

In practice, the three-antenna configuration (Fig. 13) is 
enhanced by further antenna elements so that the spac- 
ings between the antennas can be optimally adapted to 
the operating frequency range, and antenna spacings 
of a >A/2 can be used to increase the accuracy of small- 
aperture DF systems [1]. Frequently used antenna arrange¬ 
ments include the isosceles right triangle and the circular 
array (Fig. 14). 

Triangular configurations are usually restricted to frequen¬ 
cies below 30 MHz. At higher frequencies, circular arrays 
are preferred since they do the following: 
i Ensure equal radiation coupling between the antenna 
elements 

i Ensure minimum coupling with the antenna mast 
i Favor direction-independent characteristics at different 
positions due to the symmetry around the center point 


s = arccos- 


Fig. 12: Classic two-element interferometer 



Fig. 11: Principle of Doppler direction finder 



■ 


■ 
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Used in radio astronomy; block diagram (left) and antenna patterns 
(right). 
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Special consideration is to be given to avoiding ambigui¬ 
ties, which result from the fact that unambiguous mea¬ 
surement of the phase is possible only in the range of 
±180°. As already mentioned, this condition is met in the 
case of the three-element (small-aperture) interferometer 
by limiting the spacing between the elements to half the 
minimum operating wavelength. With multi-element inter¬ 
ferometers, there are the following possibilities: 
i Use of "filled" antenna arrays: Phase differences between 
neighboring elements are always < 180°; ambiguities are 
thus avoided. 

i Use of "thinned" antenna arrays: There is a phase differ¬ 
ence of > 180° between at least one pair of neighboring 
elements. Such ambiguities in antenna subarrays are 
today usually eliminated by subjecting the signals from 
all elements simultaneously to a pattern comparison by 
way of correlation > correlative interferometer. 

The basic principle of the correlative interferometer con¬ 
sists in comparing the measured phase differences with 
the phase differences obtained for a DF antenna system of 
known configuration at a known wave angle. The compari¬ 
son is performed by calculating the quadratic error or the 
correlation coefficient of the two data sets. If the compari¬ 
son is made for different azimuth values of the reference 
data set, the bearing is obtained from the data for which 
the correlation coefficient is at a maximum. 


This is demonstrated by the example of a five-element an¬ 
tenna as shown in Fig. 15. Each column of the lower data 
matrix corresponds to a wave angle a and forms a refer¬ 
ence vector. The elements of the reference vectors repre¬ 
sent the expected phase differences between the antenna 
elements for that wave angle. The upper 5 x 1 data matrix 
contains the actual measured phase differences (mea¬ 
sured vector). 

To determine the unknown wave angle (direction of inci¬ 
dence), each column of the reference matrix is correlated 
with the measured vector by multiplying and adding the 
vectors element by element. As a result, the correlation 
function K(a) is obtained, which reaches its maximum at 
the point of optimum coincidence of reference vector and 
measured vector. The angle represented by that specific 
reference vector is the wanted bearing. 

This method constitutes a special form of a beamforming 
algorithm [11], which will be discussed in greater detail in 
the following section. 


Fig. 13: Three-element interferometer 



Fig. 14: Multi-element interferometer 






Enhanced to form a multi-element interferometer. 


Fig. 15: Principle of correlative interferometer 


12 

56 

87 

I 

05 


Measured 
phase differences 




12 

56 













■ 

■ 





87 











■ 










05 






a 


Calculated phase 
differences for different 
directions of incidence 
of a plane wave 



Rohde&Schwarz Radiomonitoring&Radiolocation | Catalog 2010/2011 35 






























































Direction Finders Introduction into Theory of Direction Finding 


Direction finding using sensor array processing 
General 

The development of the classic DF methods was aimed at 
designing antenna configurations that allowed bearings to 
be determined using a circuit design as simple as possible. 
It was important to establish a simple mathematical rela¬ 
tionship between the antenna signals and the direction of 
the incident wave largely independent of frequency, polar¬ 
ization and environment. 

With the development of digital signal processing, new ap¬ 
proaches have become possible: 

i With high-speed signal-processing chips now available, 
the requirement for a simple and frequency-independent 
relationship between the antenna signals and the bearing 
no longer applies. Even highly complex mathematical 
relationships can be evaluated in a reasonably short 
period of time for determining the bearing, or handled 
quickly and economically by means of search routines, 
i Numeric methods allow the separation of several waves 
arriving from different directions even with limited 
antenna apertures (high-resolution methods, super¬ 
resolution, multiwave resolution). 

Basic design 

Fig. 16 shows a typical hardware configuration of a DSP- 
based direction finder [12]. 



Typical configuration of a DSP-based direction finder. 


The outputs of the individual antenna elements are usu¬ 
ally first taken to a network that contains the following, for 
instance: 

i Test signal inputs 

i Multiplexers if the number N of antenna outputs to 
be measured is higher than the number H of receive 
sections (tuners and A/D converters) provided in the 
direction finder 

The signals are then converted to an intermediate frequen¬ 
cy that is appropriate for the selected sampling rate of the 
A/D converters and digitized. To reduce the data volume, 
the data is digitally downconverted into the baseband. The 
complex samples x i (t) (i = 1,2, to N) of the baseband sig¬ 
nals are filtered for the desired evaluation bandwidth and 
applied to the bearing calculation section. 

Fig. 17 shows a typical implementation including a nine- 
element circular array antenna and a three-path receiver. 
The signals of the antenna elements are measured 
sequentially based on three-element subarrays. 



4 


Fig. 17: DSP-based R&S®DDF06A broadband direction finder for the 
frequency range 0.3 MHz to 3000 MHz and circular array antenna 
(R&S®ADD153SR, without cover) for the range 20 MHz to 1300 MHz 
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General DF task 

For an antenna array of N elements positioned in an 
unknown wave field, the general DF task (Fig. 18) is to 
estimate - from the measured data x 1 (t), x 2 (t), to x N (t) - the 
parameters listed below: 
i Number M of waves 
i Directions of incidence of the waves 
i Amplitudes of the waves 

Procedure 


For signal bandwidths that are small in comparison with 
the reciprocal of the signal delay between the elements 
spaced at the maximum distance (narrowband approxima¬ 
tion), the baseband signal at the output of the i-th sensor 
can be modeled in accordance with the following equa¬ 
tion: 


■ 2n, | , . 

J—|rj|cos(a 0 Pi) 

An 


x,(t)= s (t) c ,(a 0 ) e 

= s (t)a i (a 0 )+n l (t) 


+ n 


0 ) 


The task is performed in two steps: 
i Determine the relationships between the measured data 
x^t) to x N (t) and the number, directions, and amplitudes 
of the waves involved 

i Develop methods for determining the number, directions, 
and amplitudes of the main waves involved based on the 
relationships determined in the first step 

Relationship between the measured data and the 
parameters of the incoming waves 


The term s(t)= r(t)e Jp ^ 

describes the waveform and the amplitude of the signal in 
the form of the complex envelope. 

The term n, (t) 

describes the inherent noise of the sensor channel. 

The term Ci(oc 0 ) 


Data model for a single wave 

For the sake of simplicity, it is assumed that the emitter 
and the receiving antennas are located in the same plane, 
i.e. the elevation angle of the waves can be ignored. It is 
further assumed that the waves and the antenna elements 
are vertically polarized. 

For the model discussed here, it is assumed that a signal is 
radiated by an emitter with a carrier frequency of f 0 (wave¬ 
length A 0 ) and modulated with the function s(t). The wave 
is considered to be plane in the far field of the emitter and 
to arrive at an angle a 0 (Fig. 19). 


describes the characteristic of the antenna element. 

■ 271 , | , . 

j— hr os ( a o-Pi) 

The term e 0 

describes the phase displacement due to the delay 
between the reference point 0 and the position r of the 
i-th antenna element. The phase displacement solely de¬ 
pends on the position of the element (normalized to the 
wavelength A 0 ) and the direction of incidence of the wave. 


Fig. 18: General DF task 


Wave 1: s^t), a 1 Wave M: s M (t), a M 
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Measured data ^(t) x 2 (t) x 3 (t) x N (t) 

- Number of waves 
-Amplitudes 

- Directions 



Fig. 19: Characterization of incoming wave field 
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Direction Finders Introduction into Theory of Direction Finding 


In the following discussion, the antenna- and 
direction-dependent quantities are combined in the direc¬ 
tion component: 

■271, I , S 

, x j—| n| c °s(ct 0 -Pi) 

a i ( a o)= c i (a 0 )e 0 

In the interest of simplified notation and to enable straight¬ 
forward geometrical interpretation, the signals x i (t) at the 
element outputs are assumed to be components of a vec¬ 
tor x(t) in the observation space: 


x 0 ) =s 0 ) a ( a o)+ n 0 ) 


II 

-t— ' 

X 

(*.0)1 

x 2 (t) 

, n (t) = 

TOT 

n 2 (t) 


V X N OX 


V°N OX 


Characterization of antenna array 

If the wave angle a is continuously varied across the range 
of interest, the tip of the antenna vector a(a) describes a 
curve in the N-dimensional space (Fig. 20). The curve is 
referred to as an array manifold and fully characterizes the 
antenna for the parameter a, except for any loss or gain 
factors [13]. 


It is obtained by measurement or by calculation. 


Example: For an ideal crossed-loop antenna (N = 2), it 
is assumed that one element has a sine-shaped and the 
other a cosine-shaped directional pattern. The array mani¬ 
fold is expressed by the following equation: 


a (a ) = 



and thus describes a circle (Fig. 21), 


( a o) = 


^ jy~l r i l cos ( a o - Pi) ^ 

e 0 

j^T|r 2 |cos(a 0 -(3 2 ) 


2 71 

j-T|r N |cos(a 0 -p N ) 

V e J 


a(a) designates a specific direction a and is referred to as 
direction vector. The set of all direction vectors forms what 
is referred to as the array manifold, a parameter that plays 
a vital role in every aspect of array processing [13]. 


Fig. 20: Array manifold 


Element 3 



Characterizing an antenna array. 


Fig. 21: Array manifold 
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Direction Finders Introduction into Theory of Direction Finding 


Solution for single-wave model 

Assuming there is no noise, i.e. n(t) = 0, the measured 
vector x(t) differs from the direction vector a(a 0 ) only with 
respect to length. To solve the DF task, therefore, the di¬ 
rection vector has to be found, that is parallel to the mea¬ 
sured vector. The degree of parallelism can be determined 
from the direction cosine between the two vectors, which 
is proportional to the scalar product (Fig. 22). 

x a(a)= x H a(a) 

x H is the vector whose conjugate has been taken and that 
has been transposed relative to x. 


Assuming a measured vector with noise superimposed on 
it, which is the case in practice, x and a(a 0 ) are no longer 
located on a straight line. The bearing can only be estimat¬ 
ed as the best approximation between x and a(a 0 ): 


a 0 = argmax 


a 


x H a 



This approach by estimation corresponds to the beam¬ 
forming method. Like in conventional antenna arrays, the 
element signals x. are multiplied by complex weighting 
factors Wj = a, (a) and added together (Fig. 23). This 
yields a sum signal which, corresponding to the resulting 
directional pattern, depends on the direction of incidence 
a 0 of the wave and the look direction a. The asterisk (*) in 
the term above means conjugate complex. 


Fig. 22: Determining the bearing a 



The response of the output signal y to the variation of the 
weighting factors w i is used for direction finding same as 
with a classic rotating or goniometer direction finder, the 
difference being that - with numeric generation of the 
antenna pattern - the DF speed is limited only by the com 
puting power. 


Determining the bearing a 0 from the direction vector a, which is 
parallel to the measured vector x. 


Fig. 23: Beamforming by weighting the outputs of an antenna array 
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Direction Finders Introduction into Theory of Direction Finding 


If antenna arrays with largely the same elements and an 
array geometry that can be described by analytical means 
are used, the weighting factors can in most cases directly 
be calculated from the array geometry. If multiport anten¬ 
nas are used (Fig. 24), the variation of the port voltages u. 
as a function of the wave angle is as a rule determined by 
measurements. 

Since beamforming using general multiport antennas 
often does not yield a distinct directivity of the (synthetic) 
antenna pattern, the following terms are used in this case 
as well: 

i Correlation method 
i Vector matching 

Multiwave direction finding and super-resolution 
DF methods 

If unwanted waves are received in addition to the wanted 
wave in the frequency channel of interest, the conven¬ 
tional beamforming method will lead to bearing errors 
as a function of the antenna geometry. There are two ap¬ 
proaches to solve this problem: 

i If the power of the unwanted wave component is lower 
than that of the wanted wave component, the direction 
finder can be dimensioned to minimize bearing errors, in 
particular by choosing a sufficiently wide antenna aper¬ 
ture (see [1 ], chapter on multiwave direction finding), 
i If the unwanted wave component is greater than or equal 
to the wanted wave component, the unwanted waves 
must also be determined in order to eliminate them. 
When using conventional beamforming algorithms, this 
means that the secondary maxima of the DF function 
must also be evaluated. The limits are reached if either of 
the following occurs: 

■ The ratio between the primary maximum and the 
secondary maxima of the directional pattern becomes 
too small. 

■ The angle difference between the wanted and the 
unwanted wave is less than the width of the main lobe. 


By optimizing the weighting factors, the level of the 
secondary maxima can be lowered, but the width of the 
main maximum is increased at the same time. 

The super-resolution (SR) DF methods are intended as a 
solution to this problem. 

Minimum-signal direction finders are considered the 
grandfathers of the SR direction finders. In the early days 
of direction finding, bearings of co-channel signals were 
taken by alternately suppressing the waves involved by 
means of a rotating loop [1 ]. It is worth mentioning that 
signals can only be separated by audiomonitoring of the 
modulation, i.e. that correlation with acoustic patterns is 
required in order to determine the loop null. 

Beamforming methods 

Adaptive antennas are antenna arrays with beamformers 
that allow the automatic spatial suppression of unwanted 
waves [13], [14], [15]. In communications systems, this 
is done with the objective to optimize the signal-to-noise 
ratio; in direction finding, the weighting applied in order 
to suppress signals is used to determine the directions of 
incidence of the incoming waves. 

The weighting of the beamformer is selected such that 
the output power is minimized under certain secondary 
conditions. In the case of the Capon beamformer [16], the 
secondary condition for setting the weighting is defined 
such that the antenna gain remains constant for a given 
direction a. 

r 
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Direction Finders Introduction into Theory of Direction Finding 


If the incoming waves are uncorrelated, the beamformer is 
adjusted for nulls to occur in all signal directions except for 
direction a (Fig. 25). 

If the direction of an incident wave coincides with the giv¬ 
en direction a r , there is a distinct maximum in the output 
power. Fig. 26 shows an example of the angular spectrum 
of a Capon beamformer with a nine-element circular array 
(D/A = 1.4) and five uncorrelated waves. 

As with a minimum-signal direction finder, the resolution 
highly depends on the signal-to-noise ratio. 

Fig. 27 shows the same receiving scenario with noise in¬ 
creased by a factor of 10. The resolution of waves arriving 
at angles of 5° and 10° is no longer possible. 


Subspace methods 

Subspace methods are intended as a means of eliminating 
the effect of noise. This is done by splitting up the 
N-dimensional space opened up by the element outputs 
into subspaces. The common MUSIC (multiple signal clas¬ 
sification) algorithm makes use of the fact that signals 
lie perpendicular to the noise subspace. If the direction 
vectors are projected to the noise subspace, nulls that are 
independent of the noise level are obtained if signals are 
present [13], [17]. The reciprocal value is normally used as 
the DF function, so that distinct peaks occur at the direc¬ 
tions of incidence of the signals (Fig. 28 next page). 


Fig. 25: Super-resolution direction finding through nulling 



Fig. 26: DF function of Capon beamformer 


Capon beamformer 



Compared to conventional beamformer (S/N = 100); 
wave angles: 5/15/40/60/220°. 


Fig. 27: DF function of Capon beamformer 



Compared to conventional beamformer 
(S/N = 10). 
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Direction Finders Introduction into Theory of Direction Finding 


Display of bearings 

Operators of direction finders rely heavily on the display 
of DF results, which can generally be divided into two 
categories: 

i Results from a direction finder operating at a single 
frequency channel, and 
i Results from a multichannel direction finder 

Single-channel direction finders 

If a bearing is to be taken of a single frequency channel, 
the following parameters are usually displayed: 
i Bearing as a numeric value 
i Azimuth as a polar diagram 
i Elevation as a bargraph or a polar diagram 
(combined with azimuth display) 
i DF quality 
i Level 

i Bearing histogram 
i Bearings versus time (waterfall) 


Fig. 29 shows a choice of possible displays. 

In addition to the usual receiver settings such as frequency 
and bandwidth, the following parameters are set and dis¬ 
played on direction finders: 

i Averaging mode (If the signal level drops below the 
preset threshold, averaging is stopped and either 
continued or restarted the next time the threshold is 
exceeded, depending on the averaging mode.) 
i Averaging time 
i Output mode 

■ Refresh rate of display 

■ Output of results as a function of the signal threshold 
being exceeded 


Fig 28: DF function 



When using MUSIC algorithm (S/N = 10). 


Fig. 29: Bearing display with single-channel direction finding 
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Direction Finders Introduction into Theory of Direction Finding 


Multichannel direction finders 

Multichannel direction finders are implemented by means 
of digital filter banks (FFT and polyphase filters). Depend¬ 
ing on the configuration level, this enables quasi-parallel 
soil direction finding in a frequency range from a few 100 kHz 
heiften: up to a few ten MHz. Larger frequency ranges can be cov- 
einige ered by direction finding in the scan mode (Fig. 30). 

10 MHz 

With multichannel direction finders, it is essential that the 
individual events can be quickly recognized and activities 
taking place on different channels correctly assigned. 
Therefore, the following display modes are usually 
provided: 


i Bearings versus frequency 

i Bearings versus frequency and time (e.g. by displaying 
the bearings in different colors) 
i Level versus frequency (power spectrum) 
i Level versus time and frequency (e.g. by displaying the 
level values in different colors) 
i Histograms 


Fig. 30: Multichannel direction finder 
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Fig. 31: Multichannel (wideband) display 
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Direction Finders Introduction into Theory of Direction Finding 


Processing of bearings 

Position finding 

Bearings can be used in various ways to estimate the po¬ 
sition of an emitter of interest, depending on the degree 
of sophistication of the DF system and the achievable ac¬ 
curacy. Maximum accuracy is attained if several direction 
finders are employed to determine the emitter position 
by way of triangulation [1 ]. If more than two bearings are 
used for position finding, an ambiguous result is usually 
obtained (Fig. 32). 

The most probable position can be calculated in a variety 
of ways. For example, the position can be determined by 
minimizing the error squares or by maximum likelihood 
estimation. 

If a moving direction finder takes several bearings of an 
emitter, the emitter position can be determined by what is 
referred to as a running fix (Fig. 33). The emitter may move 
only slightly relative to the DF platform and must be active 
long enough to fix its position. 


Fig. 32: Position finding 


i 

i 

i 

i 



By means of triangulation. 


Fig. 33: Position finding 



By means of a running fix. 


Emitters operating in the shortwave range can be located 
by means of a single direction finder under certain condi¬ 
tions. To do this, the direction finder must be able to mea¬ 
sure azimuth and elevation, and the virtual height of the 
reflecting ionosphere layer must be known (Fig. 34). 

Emitter preclassification 

Bearings are vital characteristics of a detected signal. They 
facilitate quick segmentation, i.e. the assignment of sub¬ 
spectra to the overall spectrum of an emitter. This makes it 
possible to determine the center frequency and the band¬ 
width of an emission, which allows its automatic transfer 
to a hand-off receiver for analysis. Cluster analysis of bear¬ 
ings makes it possible to separate the signals of emitters 
operating in overlapping spectra, especially those of fre¬ 
quency hoppers. 


Fig. 35: Processing of bearings 
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DF accuracy 

DF accuracy is affected by a number of factors: 
i Wave propagation is usually disturbed by obstacles 
i Signals radiated by emitters 

■ are modulated 

■ are limited in time 

■ frequently operate at unknown carrier frequencies 

i The following are superimposed on the reception field: 

■ noise 

■ co-channel interferers 

i Noise and tolerances in the DF system 

Noise 

In addition to intermodulation distortion, interference 
caused by noise has a limiting effect on the sensitivity of a 
DF system. 

Sensitivity is understood to be the field strength at which 
the bearing fluctuation remains below a certain standard 
deviation. 


Noise may occur in the following form: 
i External noise (atmospheric, galactic, industrial noise) 
i System-inherent noise (antenna amplifiers, 

DF converters, A/D converters) 

The following basic considerations apply to system- 
inherent noise. 

Uncorrelated noise in the receive sections (Fig. 36) causes 
statistically independent variations of the measured signals 
as a function of the S/N ratio, the variations becoming no¬ 
ticeable as bearing fluctuations. 

The bearing fluctuation of general antenna arrays can be 
determined by comparing the variation of the measured 
vector x with the corresponding variation of a on the curve 
a(a). Fig. 37 demonstrates this for a two-element array. 

The "faster" the antenna vector "moves", the smaller the 
effect of the variation of the measured vector on the varia¬ 
tion of a [17]. 


Fig. 36: Model describing signal and noise 




Bearing 


The relationship between the minimum bearing fluctuation 
and the signal-to-noise ratio for a given antenna geometry 
(expressed by a(a)) is determined by the Cramer-Rao 
bound [13]. 


For a general antenna array at whose outputs K samples 
are taken, the variance o 2 of bearing fluctuation is there¬ 
fore defined by the following equation: 
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The equation simplifies if omnidirectional receiving ele¬ 
ments are used and the antenna array is of conjugate sym¬ 
metrical design: 
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Fig. 37: Relationship 



Between bearing fluctuation and signal-to-noise ratio 


For the important case of a circular array with diameter D 
and N omnidirectional receiving elements, the following 
applies at a constant elevation angle e: 
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Direction Finders Introduction into Theory of Direction Finding 


Fig. 38 shows the typical effect of the antenna diameter, 
the wavelength and the number of elements of a circular 
array on the S/N ratio that is required for a specific bearing 
fluctuation. 


If direction finding is based solely on amplitude evaluation, 
as is the case with the Watson-Watt method, the Cramer- 
Rao bound is determined by: 

2 1 f 1 1 1 

°crb_ww - 2<^SNR + 2SNR 2 , 

The effect of the S/N ratio on the bearing fluctuation is par¬ 
ticularly obvious in the case of a narrowband two-element 
interferometer. In narrowband systems, the noise voltage 
becomes approximately sine-shaped with slowly varying 
amplitude and phase so that, assuming large S/N ratios, 
the following applies to the phase variation [19]: 


9 2SNR 


Given a sufficiently long observation time, the phase varia¬ 
tion can be reduced by way of averaging. If the data used 
for averaging is uncorrelated, averaging over K samples 
improves the phase variation as follows: 



By mapping the phase variation to virtual variations of the 
positions of the DF antennas (Fig. 39), the bearing stan¬ 
dard deviation is obtained as follows: 


A 

7iD^/2SNRK 


Fig. 38: Minimum required S/N ratio 


SNR in line with the CRB 



For a bearing fluctuation (standard deviation) of 2° as a function of the 
antenna diameter as referenced to the wavelength and the number K of 
samples taken. 


For N = 2, this corresponds to the equation for g 2 rb UCA 
mentioned above, if the term 

1 

N 2 SNR 2 ' 

which can be ignored for large S/N ratios, is removed and 
s is set equal to 0. 

The above relationship clearly shows how important it is 
that the relative antenna basis D/A is as large as possible. 

Limitation of upper operating frequency because of 
ambiguity 

The required S/N ratio decreases as the frequency in¬ 
creases, because the length of the array manifold a(a) also 
increases. With a constant number N of antenna elements 
(= dimension of the observation space), however, various 
sections of the curve a(a) approach each other (Fig. 40 
next page). At the same time, the probability of measured 
vectors being positioned on such curve sections increases. 
This leads to large bearing errors [17]. 

The problem can be solved by increasing the number of 
elements N: The increase in dimension of the observation 
space thus created will accommodate the extended length 
of the array manifold. 


Fig. 39: Effect of phase noise on bearing error 
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Measurement errors 

Different gain and phase responses in the receive sections 
cause bearing errors that increase as the antenna aperture 
decreases relative to the wavelength. Fig. 41 shows this 
for a two-element interferometer. 

As already mentioned, the receive sections of most multi- 
path direction finders are calibrated for synchronization by 
means of a test generator prior to the DF operation. The 
magnitude and phase responses are measured, and the 
level and phase differences are stored. In DF operation, 
the measured values are corrected by the stored differ¬ 
ence values before the bearing is calculated. 

Special attention is to be given to the frequency response 
of the filters since synchronization must be ensured not 
only in the middle of the filter passband but also at the 
band limits. Digital filters have the decisive advantage that 
they can be implemented with absolutely identical trans¬ 
mission characteristics. 


Fig. 40: Ambiguities 



As D/A increases, various sections of the curve a(a) approach each 
other. This causes large, sudden variations of the bearing if the S/N 
ratio and the number of elements remain constant. 


Fig. 41: Effects 



Of phase synchronization tolerances on bearing error for 
different antenna apertures. 


Multi wave-related problems 

As already mentioned in the introductory section, the sim¬ 
ple case of a plane wave hardly ever occurs in practice. In 
a real environment, further waves usually have to be taken 
into account which 

i result from other emitters operating on the same 
frequency channel > incoherent interference or 
i occur in the form of secondary waves (caused by 
reflection, refraction, diffraction - see Fig. 42) > coherent 
co-channel interference (prerequisite: detours are small 
relative to the coherence lengths determined by the 
bandwidth B). 

In practice, a large number of waves is involved [20]. 

Fig. 43 shows a typical azimuth distribution of waves gen¬ 
erated by a mobile emitter in a built-up area. The direct 
wave component with amplitude 1 arrives at an angle of 
90°. 


Fig. 42: Coherent secondary waves 



Caused by reflection. 


Fig. 43: Azimuth distribution of waves 
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Radiated by an emitter in a built-up area. 
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Fig. 44 shows the resulting wavefront as a contour plot for 
the amplitude and phase [12]. 

If the majority of waves arrives from the direction of the 
emitter, the bearing error can be sufficiently reduced 
by increasing the aperture of the antenna system. This 
effect is shown by Fig. 45 for an interferometer direction 
finder. The pattern of the isophase curves in Fig. 44 shows 
that two measures are required in order to increase the 
DF accuracy: 

i The antenna aperture should be as wide as possible 
i The number of antenna elements should be significantly 
higher than with interference-free reception 

Fig. 46 shows the positive effect of a wide antenna aper¬ 
ture on the DF accuracy for a two-element interferometer. 
Arranging the antenna elements at a spacing that is large 
relative to the operating wavelength will prevent bearings 
from being obtained at points where the isophase curves 
strongly bend. 


The effect of the number of antenna elements is obvious 
when looking at the density of the isophase curves. To 
avoid 180° ambiguities in areas of high density, the spac¬ 
ing between the antenna elements should be smaller in 
these areas than with interference-free reception. This re¬ 
quires a higher number of antenna elements. 

Fig. 46 shows the effect of the antenna size and the num¬ 
ber of antenna elements on the bearing error for two cir¬ 
cular antennas in a two-wave field with an amplitude ratio 
of 0.3 of the reflected wave to the direct wave. The nine- 
element antenna delivers increasingly more accurate bear¬ 
ings up to a diameter of five wavelengths, whereas the 
bearings obtained with the five-element antenna become 
instable already from a diameter of 1.6 wavelengths. 


Fig. 44: Resulting field 



x/A 


Within a range of 4 x 4 wavelengths (the amplitudes are color- 
coded; the lines represent isophases with tt/4 spacing). 


Fig. 45: Reducing bearing error 


Incident wavs 



Of interferometer direction finder by increasing antenna aperture. 


Fig. 46: Bearing errors 



Of a five- and a nine-element circular array antenna as a function 


of the antenna diameter expressed in operating wavelengths in a 


two-wave field. 
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